---
abstractspacing: double
fontsize: 12pt
margin: 2cm
urlcolor: darkblue
linkcolor: Mahogany
citecolor: Mahogany
spacing: single
bibliography: references.bib
biblio-style: apalike
output:
  pdf_document:
    citation_package: natbib
    fig_caption: no
    number_sections: no
    keep_tex: no
    toc: no
    toc_depth: 3
    template: article-template.latex
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE, warning=FALSE, message=FALSE,cache=TRUE)
knitr::opts_chunk$set(fig.width=7, fig.height=6, out.width = '70%', fig.align = "center") 
rm(list=ls())
library(remotes)
library(kableExtra)
library(haven)
library(tidyverse)
library(ivmodel)
library(doParallel)
library(foreach)
library(estimatr)
require(AER)
library(lfe)
library(glue)
#path <- "/Users/ziwenzu/Dropbox/research/IV/IV Sensitivity/LLXZ_rep"
path <- "~/Dropbox/ProjectZ/IV Sensitivity/LLXZ_rep"
setwd(path)
knitr::opts_knit$set(root.dir = path)
#install_github("apoorvalal/ivDiag")
library(ivDiag)
# number of cores
cores <- 15
```






## @blattman2014

| Replication Summary | |
|--------|--------------|
| Unit of analysis | resident |
| Treatment | mass education campaign for dispute resolution |
| Instrument | assignment to treatment blocks |
| Outcome | serious land dispute |
| Model | Table9(8) |

```{r apsr_Blattman_etal_2014}
df <- readRDS("./rawdata/apsr_Blattman_etal_2014.rds")
df$district <- 0
for (i in 1:15) {df$district[which(df[,paste0("district",i)]==1)] <- i}
D <-"months_treated"
Y <- "fightweap_dummy"
Z <- c("block1", "block2", "block3")
controls <- c("ageover60", "age40_60", "age20_40", 
 "yrs_edu", "female", "stranger", "christian",
 "minority", "cashearn_imputedhst", "noland",
 "land_sizehst", "farm_sizehst", "lndtake_dum",
 "housetake_dum", "vsmall", "small", 
 "small2", "small3", "quartdummy", "cedulevel_bc",
 "ctownhh_log_el", "cwealthindex_bc", "cviol_experienced_bc",
 "clndtake_bc", "cviol_scale_bc", "clandconf_scale_bc",
 "cwitchcraft_scale_bc", "cpalaviol_imputed_bc",
 "cprog_ldr_beliefs_bc", "cattitudes_tribe_bc",
 "crelmarry_bc", "trainee")
cl <- "district"
FE <- "district"
weights<-NULL
(g<-ivDiag(data=df, Y=Y, D=D, Z=Z, controls=controls, FE =FE, 
  cl =cl,weights=weights, cores = cores))
```

```{r, cache = FALSE}
plot_coef(g)
```


```{r apsr_Blattman_etal_2014_sav, echo = FALSE}
save(g, file="./estimate/Blattman2014.RData")
```



